Search Results

A modified drop drag (MOD) model based on the droplet deformation and breakup (DDB) model has been developed and implemented in the KIVA II code for describing fuel spray impingement against the walls of the combustion chamber in D.I. diesel engines. The model accounts for the effects of both the drop's frontal area and its drag coefficient as a function of its distortion. In the MDD model, the stochastic rotation of drops that have impacted and bounced off the wall is considered. This makes it possible to model the highly distorted drop's frontal area variation and its effect on the drag coefficient in sprays. At the same time, the distribution of sizes of drops formed from parent drops after their breakups is modeled by using a stochastic model. A modified RNG κ - ε model is also included in the present computations. The details of fuel spray impingement against a wall were investigated by using the present models.

A numerical study of air-fuel mixing in a direct-injection spark-ignition engine was carried out. In this paper, the numerical models are described and grid generation methods to represent a realistic port-valve-chamber geometry is discussed. To model a vaporizing hollow-cone spray resulting from an automotive pressure-swirl injector, a newly developed sheet spray atomization model was used to compute the processes of disintegration of the liquid sheet and breakup of the subsequent drops. Computations were performed of a particular 4-valve pent-roof engine configuration in which the intake process and an early fuel injection scheme were considered. After an analysis of the intake-generated flow structures in this engine configuration, the spray behavior and the spatial and temporal evolution of fuel liquid and vapor phases are characterized.

A two-dimensional transient Heat Conduction in Components code (HCC) was successfully set up and extensively used to calculate the temperature field existing in real engine combustion chambers. The Saul'yev method, an explicit, unconditionally stable finite difference method, was used in the code. Consideration of the gasket between the cylinder wall and head, and the air gap between the piston and liner were included in the code. The realistic piston bowl shape was modeled with a grid transformation and piston movement was considered. The HCC code was used to calculate the wall temperature of an Isuzu ceramic engine and a Caterpillar heavy-duty diesel engine. The code was combined with the KIVA-II code in an iterative loop, in which the KIVA-II code provided the instantaneous local heat flux on the combustion chamber surfaces, and the HCC code computed the time-averaged wall temperature distribution on the surfaces.

A, three-dimensional CFD code, based on the KIVA code, is used to explore alternatives to conventional DI diesel engine designs for reducing NOx and soot emissions without sacrificing engine performance. The effects of combustion chamber design and fuel spray orientation are investigated using a new proposed GAMMA engine concept, and two new multiple injector combustion system (MICS) designs which utilize multiple injectors to increase gas motion and enhance fuel/air mixing in the combustion chamber. From these computational studies, it is found that both soot and nitrous oxide emissions can be significantly reduced without the need for more conventional emission control strategies such as EGR or ultra high injection pressure. The results suggest that CFD models can be a useful tool not only for understanding combustion and emissions production, but also for investigating new design concepts.

Advanced combustion systems that simultaneously address PM and NOx while retaining the high efficiency of modern diesel engines, are being developed around the globe. One of the most difficult problems in the area of advanced combustion technology development is the control of combustion initiation and retaining power density. During the past several years, significant progress has been accomplished in reducing emissions of NOx and PM through strategies such as LTC/HCCI/PCCI/PPCI and other advanced combustion processes; however control of ignition and improving power density has suffered to some degree - advanced combustion engines tend to be limited to the 10 bar BMEP range and under. Experimental investigations have been carried out on a light-duty DI multi-cylinder diesel automotive engine. The engine is operated in low temperature combustion (LTC) mode using 93 RON (Research Octane Number) and 74 RON fuel.

Multiple injection strategies have been revealed as an efficient means to reduce diesel engine NOx and soot emissions simultaneously, while maintaining or improving its thermal efficiency. Empirical soot models widely adopted in engine simulations have not been adequately validated to predict soot formation with multiple injections. In this work, a multiple-step phenomenological (MSP) soot model that includes particle inception, surface growth, oxidation, and particle coagulation was revised to better describe the physical processes of soot formation in diesel combustion. It was found that the revised MSP model successfully reproduces measured soot emission dependence on the start-of-injection timing, while the two-step empirical and the original MSP soot models were less accurate. The revised MSP model also predicted reasonable soot and intermediate species spatial profiles within the combustion chamber.

In spite of the efficiency of Computational Fluid Dynamics (CFD) as a design tool, numerical simulations of gaseous fuel injection have not been widely adopted because of the difficulty in modeling the complicated physical phenomena associated with high speed gas flows. In the present study, a new model for simulating transient direct injection of gaseous phase fuel, including hydrogen, into a combustion chamber using a practical computational grid was developed. The model was implemented into KIVA3V, a multi-dimensional CFD code. The new model employs several sub-models to describe the physical phenomena of high speed gas injection. The underexpanded jet issuing from the nozzle was modeled using the conditions at the Mach disk as inflow boundary conditions. The effect of turbulence is shown to lead to non-unique flow solutions.

In the Atomization regime, liquid jets breakup either within the nozzle or immediately upon entering the chamber gas and drops much smaller than the jet diameter are formed. The mechanism of Atomization, which is presently unknown, was investigated by the simultaneous use of two photographic techniques. The initial transient was observed with a 106 frames/s camera and the steady state by a technique similar to spark photography. The experiment range was: liquid pressure 500 to 2500 psia; five mixtures of water and glycerol to vary the liquid viscosity; air, nitrogen, helium, and xenon at up to 600 psia as chamber gases to separate gas pressure from gas density effects; and 14 nozzle designs. Not changed were the temperature (room value), the nozzle diameter (340 μ), and the surface tension (70 dyne/cm).

The effect of auxiliary gas injection on diesel engine combustion and emissions was studied using KIVA, a multidimensional computational fluid dynamics code. Auxiliary gas injection (AGI) is the injection of gas, in addition to the fuel injection, directly into the combustion chamber of a diesel engine. The objective of AGI is to influence the diesel combustion via mixing to reduce the emissions of pollutants (soot and NOx). In this study, the accuracy of modeling high-speed gas jets on very coarse computational grids was addressed. KIVA was found to inaccurately resolve jet flows near walls. The cause of this inaccuracy was traced to the RNG k - ϵ turbulence model with the law-of-the-wall boundary condition used by KIVA. By prescribing the length scale near the nozzle exit, excellent agreement between computed and theoretical jet penetration was obtained for a transient jet into a quiescent chamber at various operating conditions.

A phenomenological nozzle flow model has been developed and implemented in both the FIRE and KIVA-II codes to simulate the effects of the nozzle geometry on fuel injection and spray processes. The model takes account of the nozzle passage inlet configuration, flow losses and cavitation, the injection pressure and combustion chamber conditions and provides initial conditions for multidimensional spray modeling. The discharge coefficient of the injector, the effective injection velocity and the initial drop or injected liquid ‘blob’ sizes are calculated dynamically during the entire injection event. The model was coupled with the wave breakup model to simulate experiments of non-vaporizing sprays under diesel conditions. Good agreement was obtained in liquid penetration, spray angle and drop size (Sauter Mean Diameter). The integrated model was also used to model combustion in a Cummins single-cylinder optical engine with good agreement.

An electronically controlled Caterpillar single-cylinder oil test engine (SCOTE) was used to study diesel combustion. The SCOTE retains the port, combustion chamber, and injection geometry of the production six cylinder, 373 kW (500 hp) 3406E heavy-duty truck engine. The engine was equipped with an electronic unit injector and an electronically controlled common rail injector that is capable of multiple injections. An emissions investigation was carried out using a six-mode cycle simulation of the EPA Federal Transient Test Procedure. The results show that the SCOTE meets current EPA mandated emissions levels, despite the higher internal friction imposed by the single-cylinder configuration. NOx versus particulate trade-off curves were generated over a range of injection timings for each mode and results of heat release calculations were examined, giving insight into combustion phenomena in current “state of the art” heavy-duty diesel engines.

A study was performed that takes into account both turbulence and chemical kinetic effects in the numerical simulation of diesel engine combustion in order to better understand the importance of their respective roles at changing operating conditions. An approach was developed which combines the simplicity and low computational and storage requests of the laminar-and-turbulent characteristic-time model with a detailed combustion chemistry model based on well-known simplified mechanisms. Assuming appropriate simplifications such as steady state or equilibrium for most of the radicals and intermediate species, the kinetics of hydrocarbons can be described by means of three overall steps. This approach was integrated in the KIVA-II code. The concept was validated and applied to a single-cylinder, heavy-duty engine. The simulation covers a wide range of operating conditions.

Combustion simulations utilizing the characteristic time combustion model have been performed for four DI diesel engines ranging in size from heavy-duty to large-bore designs. It has been found that the pre-factor to the turbulent characteristic time acts as a scaling parameter between the engines. This phenomenon is explained in terms of the non-equilibrium behavior of the turbulent time and length scales, as is encountered in the rapidly distorting, spray-induced flows of DI diesel engines. In fact, the equilibrium assumption between turbulence production and dissipation, which forms the basis for the employed k-ε-type turbulence models, does not hold in these situations. For such flows, the real turbulent dissipation time scale is locally proportional to the turbulent characteristic time scale which is determined by a typical eddy turnover time.

The effects of numerical methodology in defining the initial conditions and simulating the compression stroke in D.I. diesel engine CFD computations are studied. Lumped and pointwise approaches were adopted in assigning the initial conditions at IVC. The lumped approach was coupled with a two-dimensional calculation of the compression stroke. The pointwise methodology was based on the results of an unsteady calculation of the intake stroke performed by using the STAR-CD code in the realistic engine and port geometry. Full engine and 60 deg. sector meshes were used in the compression stroke calculations in order to check the accuracy of the commonly applied axi-symmetric fluid dynamics assumption. Analysis of the evolution of the main fluid dynamics parameters revealed that local conditions at the time of injection strongly depend on the numerical procedure adopted.

An experimental study of luminous combustion in a modern diesel engine was performed to investigate the effect of injection parameters on NOX and soot formation via flame temperature and soot KL factor measurements. The two-color technique was applied to 2-D soot luminosity images and area-averaged soot radiation signals to obtain spatially and temporally resolved flame temperature and soot KL factor. The imaging system used for this study was based on a wide-angle endoscope that was mounted in the cylinder head and allowed different views of the combustion chamber. The experiments were carried out on a single-cylinder 2.4 liter D.I. diesel engine equipped with an electronically controlled common-rail injection system. Operating conditions were 1600 rpm and 75% load. The two-color results confirm that retarding the injection timing causes lower flame temperatures and NOX emissions but increased soot formation, independent of injection strategy.

Intake flow simulations were carried out for a prototype DISI engine using the standard k-ε model and the RNG k-ε model. The results were compared with PTV (transient water analog) measurements. The study was focused on low load operations with engine speed at 400 rev/min. Two cases were studied, a single intake case in which one intake port was blocked and a dual intake port case. In the computations, the results show that the standard k-ε model tends to produce higher turbulence levels when turbulence is generated and decays faster when turbulence dissipates. Different turbulence models predict almost the same flow structures. However, the effects of the turbulence model on the predicted tumble and swirl ratios are significant. The TKE distributions at BDC predicted by the two models are also different. The standard k-ε model seems to be more diffusive. Good agreements with PTV data were obtained in the single valve case with the RNG k-ε model.

Fuel preparation and stratified combustion were studied for a conceptual gasoline Direct-Injection Spark-Ignition (GDI or DISI) engine by computer simulations. The primary interest was on the effects of different injector orientations and the effects of tumble ratio for late injection cases at a partial load operating condition. A modified KIVA-3V code that includes improved spray breakup and wall impingement and combustion models was used. A new ignition kernel model, called DPIK, was developed to describe the early flame growth process. The model uses Lagrangian marker particles to describe the flame positions. The computational results reveal that spray wall impingement is important and the fuel distribution is controlled by the spray momentum and the combustion chamber shape. The injector orientation significantly influences the fuel stratification pattern, which results in different combustion characteristics.

To model sprays from pressure-swirl atomizers, the connection between the injector and the downstream spray must be considered. A new model for pressure-swirl atomizers is presented which assumes little knowledge of the internal details of the injector, but instead uses available observations of external spray characteristics. First, a correlation for the exit velocity at the injector exit is used to define the liquid film thickness. Next, the film must be modeled as it becomes a thin, liquid sheet and breaks up, forming ligaments and droplets. A linearized instability analysis of the breakup of a viscous, liquid sheet is used as part of the spray boundary condition. The spray angle is estimated from spray photographs and patternator data. A mass averaged spray angle is calculated from the patternator data and used in some of the calculations.

A correction for the turbulence dissipation, based on non-equilibrium turbulence considerations from rapid distortion theory, has been derived and implemented in combination with the RNG k - ε model in a KIVA-based code. This model correction has been tested and compared with the standard RNG k - ε model for the compression and the combustion phase of two heavy duty DI diesel engines. The turbulence behavior in the compression phase shows clear improvements over the standard RNG k - ε model computations. In particular, the macro length scale is consistent with the corresponding time scale and with the turbulent kinetic energy over the entire compression phase. The combustion computations have been performed with the characteristic time combustion model. With this dissipation correction no additional adjustments of the turbulent characteristic time model constant were necessary in order to match experimental cylinder pressures and heat release rates of the two engines.

In this work computational and experimental approaches are combined to characterize in-cylinder flow structures and local flow field properties during operation of the Sandia 1.9L light-duty optical Diesel engine. A full computational model of the single-cylinder research engine was used that considers the complete intake and exhaust runners and plenums, as well as the adjustable throttling devices used in the experiments to obtain different swirl ratios. The in-cylinder flow predictions were validated against an extensive set of planar PIV measurements at different vertical locations in the combustion chamber for different swirl ratio configurations. Principal Component Analysis was used to characterize precession, tilting and eccentricity, and regional averages of the in-cylinder turbulence properties in the squish region and the piston bowl.